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■ Abstract 

CO , We present exact solutions of a class of the nonlinear models which describe the parametric 

' conversion of photons. Hamiltonians of these models are related to the classes of finite orthog- 

(-H I onal polynomials. The spectra and exact expressions for eigenvectors of these Hamiltonians are 

^ |. obtained. 
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> 

H ■ 1 Introduction 

a 

In nonlinear optical models the influence of a medium on electromagnetic field (^E , is described 

by the material source-free Maxwell equations, where in general a functional dependence P = P [-E] 

of the polarization P on the electric field E is assumed (see e.g. jP-L| . |H-0-T| . |B-Cp . This 
dependence describes complicated microstructure of the medium and the nonlinearity of the matter - 
field interactions. Assuming the classical description of the medium, the field can be quantized, and 
the energy operator is obtained. In the case of the two-mode field the operator is given by 

H = Ho + e"'"" * H/ e'^" \ (1.1) 

where 

Ho = woaoao + tJia;jai (1.2) 
describes the free field, and the term 

oo 

H/= {c^Mmn{^f{^lf^S.^ + h.c), (1.3) 
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where aioio = ceoioi = 0; is the interaction Hamiltonian responsible for the Ught-matter interactions. 
The annihilation and creation operators of two modes ao,ai,a|5,a* fulfill the Heisenberg canonical 
commutation relations. Using the boson-number ordering, see |0-H-Tj . H/ can be expressed in the 
form: 



H/=| E /feofci(aSao,aJai)ag«a^+/i.c. I + (1.4) 

oo 

9kok^ (a^ao, a*ai)ago (a*)''' + h. 

,kQ,ki—0 

where fkoki{x,y) and gkoki(x,y) are functions of two arguments x, y. These functions are defined by 
the constants akimn and are responsible for the light-matter interaction via the functional dependence 
P = P[E]. If the sum in H1.3|l is finite then fkoki{x,y) and gkoki{x,y) are polynomials. This case 
has been investigated during the last decade by many authors, see e.g. |(t-K-()| . |,)-D| . They 
used approximate or semiclassical methods. 

Let us give the interpretation of the particular terms of the Hamiltonian l|1.4(l . The term 
gkoki (agao, a*ai)aQ'' (a^;) ~^ describes the process of simultaneous absorption of fco photons in mode 
and emission of fci photons in mode 1. The probability of this process depends on g^^fej^ (a^ao, a|fai), 
i.e. it depends on the intensity of light in the medium . So the function gkoki is a generalization of 
the coupling constant for the conversion. Such a process is called the intensity dependent |,T-D| (or 
parametric |Pe-Lu| ) conversion of feg photons in the mode into fci photons in mode 1. The hermitian 
conjugate term (^gk^kii^o^OT ^i^i)^o° (^i)'^^ )* describes the parametric conversion of fci photons in 
mode 1 into fco photons in mode with coupling given by the operator g^-^^j^ i^o^ + ^Oi ajai — fci). By 
analogy the term /^^fc^ (apao, a*ai)aQ°a^^ corresponds to the process of absorption by the medium of 
the cluster consisting of the fco photons in the mode and fci photons in the mode 1. The hermitian 
conjugate term describes the emission of the same cluster by the medium. 

In this paper we study the intensity dependent conversion of the fixed number fco of photons in 
mode into a fixed number fci of photons in mode 1 and vice versa. To simplify the notation let 
Hq := goQ + 5oo ^^'^ 9 gkoki', the interaction Hamiltonian for such process takes the form 

H/ = /io(a*ao,aJai) +g(a*ao,aJai)aS« (aj)''^ + (a*)*^" aj^^g (a^ao, ajai) (1.5) 

with the complex-valued function 

g (ajao, ajai) = e^^Kaca^a,) |^ (a^ao, ajai) |. (1.6) 

This paper is a continuation of the program initiated in j()-H-T| . |H-()-T| . where the theory 
of orthogonal polynomials was applied to the solution of the eigenproblem for the Hamiltonians of 
the type (|1.5|l for general functions ho and g. The solution of the eigenproblem of the interaction 
Hamiltonian H/ is equivalent to the problem of the integrability of the system under consideration; 
it follows from the fact that the solution of Schrodinger equation for the Hamiltonian (|1.1() is 

In Section 2 a family of operators commuting with the Hamiltonians H1.2|l and (|1.8(l is found, i.e. 
the integrals of motion of the system. It is shown that the Hilbert space can be split into finite- 
dimensional subspaces invariant under the action of Ho and H/. Introducing the Fock basis in the 
subspaces, the matrix form of the interaction Hamiltonian is the Jacobi matrix. In Section 3 it is 
shown that the eigenproblem of the considered Hamiltonian is equivalent to the moment problem of the 
theory of finite orthogonal polynomials. Using this equivalence we express the spectral decomposition 
of the Hamiltonian in terms of the orthogonal polynomials. Section 4 is devoted to the construction 
of the families of Hamiltonians related to the fixed systems of orthonormal polynomials. In Section 5 
the families of integrable Hamiltonians are presented. We present examples of simplest Hamiltonians 
related to the well known families of finite orthogonal polynomials. The spectra and eigenvectors of 
these Hamiltonians are presented. 
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2 Conversion of fixed numbers of photons 

Let us show that the two mode Fock space TL can be split into finite dimensional subspaces, invariant 
with respect to the interaction Hamiltonian (|1.5|) and the free Hamiltonian \i.2\ . and therefore also 
invariant with respect to the Hamiltonian First we observe that the following operators 

K := fcia^ao + /coata^ (2.1) 



R.:=^ + E ZJi:^' ^ (2-2) 



where k = 0, 1, commute with Hq and H/, i.e. they are integrals of motion of the system, fco and fci 
are defined by (|1.5|l . The elements of Fock basis of Ti 

|no,ni) = ^i==(a5)"''(at)"M0,0), no, m =0,1,... (2.3) 

are eigenvectors of the operators (|2.1|l . H2.2|l : 

K |no,ni) = (fcino + fconi) |7io,ni), (2.4) 
|f^o,f^l) = f« |"o,"-i) , (2.5) 

where the eigenvalues r„ are equal to the remainder of the division of by {r^ = nf^{'n\oA /ck)), 
see | G-K-0| . So one has the orthogonal decomposition 

W = 0Wm (2-6) 

of TL into finite-dimensional Hilbert subspaces labelled by multi-indices /i := (ro,ri,iV) e J := 
{0, 1, . . . , fco ^ 1} X {0, 1, . • . , fci — 1} X (N U {0}) and spanned by the vectors 

|n)^ := |ro + fcon, ri + fci(iV-n)) n = 0, l,...,iV (2.7) 

of the Fock basis. The Hilbert subspace is a common eigenspace of the operators K, Rq and Ri 
with dimension 

dimTi:^ = iV+ 1. (2.8) 
Moreover, + 1 can be obtained as the only eigenvalue of the dimension operator 

D = -^K - -^Ro - -^Ri + 1 (2.9) 

kqKi ko ki 

on Tip. According to the above let us introduce the operators 

Ao := ^a^ao, (2.10) 

A:=g(a5ao,atai)ag°(a|)'^ (2.11) 

and also replace the operators a^aa and sl\sli by Aq and K. 
The operators Aq, A, A* satisfy the relations 

[Ao,A] = -A, [Ao,A*]=A*, 

A*A = g(Ao-l,K), (2.12) 
AA* =g(Ao,K), 

where the function Q (Aq, K) is determined by g: 



e(Ao,K) 



(fcoAo + l)(fcoAo + 2) . . . (fcoAo + fco) x 



g (^fcoAo, -^K - fciAo 
i-K - fciAo") ("^K - fciAo - iV . . (^K - fciAo - fci + 1 ) . (2.13) 
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It can be shown that 
ko 



Ao\n)^^ {^+n]\n)^, A|n)^ = 6„_i,^ |n - 1) A* = 6„,^ jri + 1) (2.14) 



where 

h 



V-o / Wt/ I — +n, fciro + fcori +fcofciAf I, (2-15) 



i.e. Ao is diagonal, A and A* are weighted shift operators and (j) denotes the phase factor of function 
g expressed in the new variables n and N. Let us notice that from ()2.13f) it follows that ^ (— 1, N) = 
Q [N, N) = what makes ^^TWji consistent. 

The interaction Hamiltonian H1.5|l now takes the form 

H/ = /i(Ao,K)+A + A*, (2.16) 

where h is uniquely determined by ho and the free Hamiltonian Hp is given by 

Ho = {ujoko - c^ifci)Ao + "^K. (2.17) 

Thus the Hamiltonian describing our system belongs to the operator algebra generated by the oper- 
ators K, Aq, a, a* (K commutes with the others). The subspaces are invariant subspaces of 
the operators Hq, H/ and therefore of H. Moreover, the action of these operators on elements of the 
Fock basis is 

Ho \n)^ = (^uoTq +cjiri + ujikiN + (cjo^o - u;iki)n^ \n)^ (2-18) 

and 

H/ = bn-i,fj. \n - 1)^ + a„,^ \n)^ + bn,p. \n + 1)^, (2-19) 
where bn,^, is given by H2.15|l . and 

a„,^ = + n, kiro + kon + kokiN^ . (2.20) 

It follows from l|2.19|l that the matrix form of the operator := H/ in the basis { 1"-)^ }n-o 
Ti^i is the Jacobi matrix (three diagonal and hermitian matrix). This fact allows us to use the theory 
of orthogonal polynomials to solve the eigenproblem of |0-H-T| . |H-0-T| . 



3 Integration and finite polynomials 

The three-term formula H2.19I) suggests that the eigenproblem of is strictly connected with the 
theory of finite orthogonal polynomials. One can apply this theory under the additional assumption 
that Hp has + 1 different eigenvalues {i?;,p};_p with the corresponding eigenvectors { |i?z,p) };_o' 
i.e. 

n^\Ei,^) = Ei^,, \Ei^^). (3.1) 
Then the spectral decomposition of the interaction Hamiltonian is 



If we decompose \Ei^^) in the Fock basis {|?^)p}^^=o "^m 

N 

\Ei,f.)^Y.PniEL^^) I") A" (3-3) 



ji=0 
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then from l|2.19|l . (|3.1|l and H3.3|l it follows that the coefficients P^^Ei^^) satisfy the three term identity 



Ei,^. Pr^XEl,^.) = PiU{Ei,^,)+an^^ Pl^iEu^.) + bn.^ n,l = 0,1, N (3.4) 

which can be considered as the recurrence relations for PI^{Ei^^) with the initial condition Pq{Ei^^)= 1. 
Thus PI^{Ei^fj,) is a polynomial of degree n in the variable i^;,^. Since is hermitian, the set 
Ai)}i^o forms an orthogonal basis in 7i^. The orthogonality relations for eigenvectors \Ei^fj_) imply 
the orthonormality relation in the set of polynomials {PniEi^fj,)}'^^^: 

N 1 

^ P„^(£;,,^) Pl^^iEi^^)——— = Smn ; (3.5) 
^ {Ei,f,\Ei^^,) 

this allows to invert the formula (|3.3|l : 

In such a way we obtain a finite system of orthonormal polynomials <^PI^j^_^, with respect to 
the weight function — dependent on the discrete variables |£';^^}^p. 

Since the interaction Hamiltonian H/ preserves the decomposition ()2.6|l , the unitary one-parameter 
group of Schrodinger evolution in the interaction picture e~^^' ' (see (|1.7ll ') preserves the decomposition 
of the identity 

N 

Mi'^l- (3.7) 

Thus, using the spectral decomposition H3.2(l and the orthogonal relations (|3.5|l the operator e~'^^ * 
can be expressed in the form 

,vL,t ^yy ^-^E,,,t \Ei,^.) {Ei,^\ ^ ^ ^ 

Its matrix elements 

Am\e-'^' * |n)^ = Sr,s -^.o^o^ng. E ^™(^'^a<) ^'^(^'^A') (e \E V ^^'^^ 

where — {ro,ri,N) and = (go, 91, 5"), are expressed in terms of the orthogonal polynomials. 
From l|1.7() and H3.9|l the time evolution of the expectation value of any quantum observable X in a 
normalized state £ H becomes 

oo N S 

{X{t))^ ^ (V.|e^H^*e'"°*Xe-^««*e-^H^* IV^) = E E E ('/'l"^)^' ^(HXk). .(n|V')x 

fj,,u^J m,r,/— /c,s,n— 



—it I (.jo(go-n))+wi(9i— J'i)+[»jifci (S— A'') + ((.jofeo -1^1 fci )(«—'") 
X e ^ ^ X 



g -zt(iSfc,,-i5i,^) 



^ IF N/F IF N ^™(^^;,m) PnEk,.)PnEi,^) P^iEk,.). (3.10) 
Similarly matrix elements of X(t) are equal 

^(m|X(t) = 2^ 2^ ^(r|X|s)^ e V 

i,r=0 fc,s=0 

-«(£;.,.„-£;,„) 



e 



^ /F IF \/F IF \ ^™(^'.m) Pi'(i?fc,.)TO,M) ^»,.)- (3.11) 
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4 Hamiltonians related to the same families of fixed orthonor- 
mal polynomials 

It is natural that different Hamiltonians after reduction can lead to the same family of orthogonal 
polynomials. In this section we solve the inverse problem, i.e. how to construct different Hamiltonians 
which are related to the same family of orthogonal polynomials. 

Let us consider the interaction Hamiltonian in the case when kg = fci = 1 

H/ = /lo (aoao,a;fai) + .g (apao, a^fai) aoaj + a^aig (aoao, a^ai) . (4.1) 

Here the multi-index /i describing the decomposition H2.6|l is a single index fi = N E NU{0}. Moreover, 
all components of that decomposition have different dimensions. Let us assume that the eigenproblem 
is solved by a family of orthonormal finite polynomials {Pn {Ei,N)}n=o- We will call H4.1II the initial 
Hamiltonian. 

For any fixed fco, /ci e N we construct the Hamiltonian 

Hj = /^o(a*ao,aJai) +g(a*ao,aJai)ag° (a^)^^ + (a5)'=° aM (a^ao, ajai) (4.2) 

in such a way that the eigenproblem for any H^t, where fi = (ro, 7'i, A^) (z J — {0, 1, . . . , fco — 1} x 
{0, 1, . . . , fci - 1} X (N U {0}), is solved by the initial family {P^ {Ei.n)}^^q, i.e. 

Wfi = (ro, ri, TV) e J Pi: = P.^ and Ei^^ = E^n, 7i, ^ = 0, 1, . . . , TV. (4.3) 

In other words, the solution of eigenproblem for does not depend on ro and ri. 
Let us introduce the following operator 

I¥fcofe,(a*ao,a*ai) ^ (aga^ - Ro + fco)(atai - Ri)x 

X ^ =. (4.4) 

v/fcofci(a5ao + l)(a5ao + 2) . . . (agap + ka)a.lai{a.laii - 1) . . . (a^fa;^ - fci + 1) 

In particular VFii(aQaQ, a^a^) = 1. It is easy to check that 

Wu,uA^2.^,<^^)4"i<f' l">p ^^/n{N-n+l)\n-l)^. (4.5) 
Since additionally we have that 

^(a5ao-Ro)|n)^ = n|,^)^ (4.6) 

^(a|ai-Ri)|n)^ = (iV-n)|n)^, (4.7) 

then putting 



and 




we obtain the Hamiltonian H4.2(l which satisfies the conditions H4.3|l . 

In the next section we present some examples of initial Hamiltonians which are related to the 
families of the finite discrete orthonormal polynomials, which are best known in the literature. 

5 Quantum systems related to some classes of finite orthonor- 
mal polynomials 

In this section we present a list of Hamiltonians for which the spectral decompositions can be expressed 
by some selected families of finite orthogonal polynomials. The results of the previous section allows 
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to restrict our list to the cases when ko = ki — 1. The notation used in this section is the same as in 
|K-Sj . The hypergeometric series are denoted by 



p I Cll , • • • , Ctr 

bi,...,bs 



E 



(ai, . . . ,ar)n z'^ 



^0 ibi^---^bs)n n\' 



(5.1) 



where (ai, . . . , ar)n ■= {ai)n • • • (ar)n, ia)n ■= o,{a + l)(a + 2) • • • (a + n — 1) for ti = 1,2,... and 
(a)o ■— 1. The basic hypergeometric series is defined by 



fli, . . . , Or 

6i, ... ,6s 



q; z 



E 



. . .,bs;q)n 



(5.2) 



for < g < 1, where (ai, . . . ,ar;q)n := (ai;g)n • • • {ar]q)n and (a;g)„ := (I - a){l- aq){l- aq^) ■■■{!- 
aq'^~^), for ti = 1, 2, . . . and (a; (7)0 = 1. 

5.1 Integrable systems related to the Krawtchouk polynomials 

The Krawtchouk polynomials arise for a system described by the Hamiltonian 



Hi =p a*ai + (1 - p) a^ao + vW-p)( aoaj + a^a^), 
where < p < 1. The spectrum of H/ is 

a(H/) =N U{0}, 



(5.3) 



(5.4) 



which follows from the fact that the eigenvalues of the reduced Hamiltonian Hjv are Ei^n — I for each 
N. Thus the eigenspaces Ti.\ I E cr(H/) are the infinite-dimensional Hilbert subspaces 



= span||i;i,w) : N ^ I + n, n = 0, 1,...|, 



where 



N 



\Ei,n) ^Y.^"^^' 

^n;p,N) \n, N -n) 



(5.5) 



(5.6) 



and 



(-l)"n! \l-p 



2F1 



—n, —I 
-N 



(5.7) 



are the Krawtchouk polynomials. For each fixed N, due to 13.511 . the finite family < KniEi n',P, N) > 
forms an orthonormal system with respect to the weight function 



1 



N 
I 



p'(l-p) 



N-l 



We can summarize that the spectral decomposition of the interaction Hamiltonian is 

00 N 



H. = EE^ 



N=0 1=0 



(5.8) 



(5.9) 



The Hamiltonian (|5.3() is quadratic in annihilation and creation operators and therefore the system 
is also integrable via Heisenberg equations. 
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5.2 Integrable systems related to the Dual Hahn polynomials 

The Dual Hahn polynomials are related to the system given by the Hamiltonian 

H/ = aoao(aiai + (5+ 1) + (agao + 7 + l)atai + (5.10) 



(a*ai +(5)(a*ao + 7+ l)aoa:j + a5ai^(a*ai + (5)(a*ao + 7 + 1) 

where 7, i5 > — 1. The spectrum of this Hamiltonian is 

cr(H/) = + 7 + 5 + 1) : / = 0, 1, . . . } (5.11) 

and the eigenvalues of the reduced Hamiltonian H^r are Ei,n = Z(/+7 + (5 + l) and do not depend on N. 
For each eigenvalue l{l + 'f + 6 + 1) the corresponding eigenspace of H/ is the infinite-dimensional 
Hilbert space 

= span||£;i,jv) : N ^ I + n,n = 0,1, .. (5.12) 



where 



N 



\Ei^n) = J2 Rn{Ei,N\l, S, N) \n, N~n) (5.13) 



n=0 



and 



RAE,^; 7, s, N) = ,/ (7+") (^t^r") ( ^; ! + ' + ' 



(5.14) 



are the Dual Hahn polynomials. The finite family i. Rn(Ei n;j,6, N)\ , due to 13.5II . forms an 
orthonormal system with respect to the weight function 

1 _{2l + j + S+l){j + l)i{-iy{-N)im ^^^^^ 



{Ei,n\Ei,n) (Z+7 + <5+1)jv+i(-^ + 1); 

We can summarize that the spectral decomposition of the interaction Hamiltonian is 

00 N 



(5,16) 

5.3 Integrable systems related to the discrete Chebyshev polynomials 

The Hamiltonian 

rr ^ (2agao + a*ai + l)agao (a*ao + l)a*ai 
' 2(2a*ao + l) ^ 2(2a5ao + 1) ^ ^ ^ 

1 / (2agao + a^ai + 2)(aaao +T) ^ a* + a*a - / (^^Sao + a^ai + 2)(aSao +T) " 
2Y (2a*ao + l)(2aSao + 3) '''>''^2^ (2aSao + l)(2aSao + 3) 

describes the system for which the solution of the eigenproblem is given by discrete Chebyshev poly- 
nomials. The spectrum of (|5.17|1 is 

ct(H/)=NU{0} (5.18) 

and similarly as in subsection 5.1, the infinite-dimensional eigenspaces of H/ are spanned by the 
vectors 

N 

=^T„(£;,,jv;iV)|n,Ar-n>, N = I + n, n ^ 0,1, . . . (5.19) 

n=0 
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with the eigenvalues Ei^n = h Z = 0, 1, ... of the reduced Hamiltonian Hjv- The coefficients 

T(F -m- / (2" + lj(--V)«A^! ^ f ~n,n+l,-l, \ 



are the discrete Chebyshev polynomials. The weight function for the family |t„(£';^jv; ■^)| is, as 

I- ' J n=0 

for classical Chebyshev polynomials, constant 

= 1. (5.21) 



{Ei^n\Ei^n) 

We can summarize that the spectral decomposition of the interaction Hamiltonian is 

oo N 



H/= 5^5^^ \Ei,n) {Ei,n\. (5.22) 

JV=0 1=0 

5.4 Integrable systems related to the Hahn polynomials 

The Hahn polynomials arise for the system described by the Hamiltonian 

jj ^ agao(2agao + a*ai+a + /?+l)(aSao+/?) 
^ (2a5ao + a + /3)(2aSao + a + /3+l) ^ ' ' 

(aj^ao + a + l)(a;5ao + a + /3 + l)a5;ai _^ 



(2a*ao + a + l3 + l)(2a*ao + a + /? + 2) 



/ (2agao + a*ai + a + f3 + 2)(agao +(3 + l)(a*ao + a + l)(a*&Q + a + (3 + I) 
V (2aSao + a + /3+l)(2a*ao + a + /3 + 2)2(2a5ao + « + /? + 3) 



+ agaii 



' (2a5ao + ajai +« + /? + 2)(a5ao +(3+ l)(a5ao + a + l)(a2ao + a + p + l) 



(2a5ao + a + f3+ l)(2a5ao + a + p + 2)2(2aSao + a + /3 + 3) 
where a, /? > — 1. The spectrum of H/ is 

C7(H/) = N U {0}. (5.24) 
and the infinite dimensional eigenspaces are spanned by the vectors 

N 

\Ei,N)=Y,Qn{Ei,N;a,f3,N)\n,N-n), N = I + n, n = 0,1, . . . , (5.25) 
where Ei^n = I and 

/ (Or, J- ^ -L R M ^^^^ -L ^\ (—n\ m / _^ ^ _l « _l i \ 

1 (5.26) 



n tjr ^ « ^r^ (2n + a + /3 + l)(a + l)„(-iV)„ iV! ^ f -n,n + a + (3 + 1, -I, 

V (-l)"(n + a + /3 + l)jv+i(/3 + l)„n! V a + 1,-N 



are the Hahn polynomials. The weight function for the family \Qn(Ei N;a,p,N)\ is 

I. ' J n=0 

1 / a + l \ ( p + N-l 



{Eln\Ei,n) \ I J\ N-l 
We can summarize that the spectral decomposition of the interaction Hamiltonian is 

oo N 



(5.27) 



\Ei,n} {El 
{Ei^n\Ei^ 

Let us note that putting a = /? = we obtain the discrete Chebyshev polynomials 



^ot^o {Ei,n\Ei,n) 
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5.5 Integrable systems related to the Dual g— Hahn polynomials 

For any fixed < g < 1 and < 7, 5 < such that < -yd < the Hamiltonian 

Hj = l+jSq--fq{l-q^>''){S-q-^>'-'^)-{l-q-^>'){l--fq^>''+'^) 



(5.29) 



+ 



(a^ao + l)a;fai 



-aoai 



(a^ao + l)a^ai 
has the spectrum 

a(H^) = {g-'+7V+':/ = 0,l,---} 
and the eigenspace corresponding to eigenvalue + 7(5^'+^ e o-(H7) is 

= span||£;i,Ar) : TV = ^ + n, n = 0, 1, . . . | 

with 

-E;,jv = g-'+7V+', 

and 



|i^(,]v) = ^ Rn{Ei,N-n, S, N\q) \n, N - n) 



n=0 



and 



' (^g;g)iv(7g,9-^;g)n(7g)^ 
(7^g2;g)Ar(g,5-lg--^;g)„(7(5g)" ^ ^ 



q ",q ',7(5g'+^ 

— AT 



(5.30) 



(5.31) 



(5.32) 



(5.33) 



g;g (5.34) 



which are the Dual q Hahn polynomials. For each fixed N the family ^-E(,jv; 7, -^|g) | forms 
an orthonormal system with respect to the weight function 



1 



(79,7^^9, g-^;g),(l-7<5g^'+^) 



{Ei,n\Ei,n} {q,j6q^+^,6q;q)i (1 - 7(5g) (-7g)' 
We can summarize that the spectral decomposition of the interaction Hamiltonian is 

00 N 



\^i,n\-i^i,n) 



(5.35) 



(5.36) 



N=0 1=0 



5.6 Integrable systems related to the AfRne g— Krawtchouk polynomials 

For any fixed < g < 1 and < p < g^^ the Hamiltonian 

1 - [(1 - q-'''^'^'){l - pq''°'^+^) - pq-''^^ (1 - q''>'') 



H7 



(5.37) 



+ 



_pq-a.la.^ + l^l _ ga2ao + l)(-X - g-*^*i)(l - pq^o^o + l^ 



(agao + l)aiai 



* I 

aoa^ + 



+ aoan 



has the spectrum 



(agao + l)a^ai 
<7(H,) = {g-' :Z = 0,1,...}. 



(5.38) 
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Eigenspaces H , q € cr(H/) are 



= sj9an||£';,jv) : N = I + n, n = 0,1, . . 



for Ei^N = q 



N 



\Ei,n) = Y.^n"iEi,N;p,N;q)\n,N -n) 



ra=0 



and 



V (9;9)n(g;9)jv-n V pI'I 



(5.39) 



(5.40) 



(5.41) 



which arc the affinc (/-Krawtchouk polynomials. The family ^K^^-f (Ei^n;p, N; q)^ is the or- 
thonormal system with respect to the weight function 



(pq-.q)i{q;q)N , .-i 
7 — w — ^ — (Pi) ■ 



{Ei^n\Ei^n) {q;q)i(q;q)N-l 
We can summarize that the spectral decomposition of the interaction Hamiltonian is 

oo N 



JV=0 1=0 



-i \Ei,n) {Ei,n\ 
{Ei^n\Ei^n) 



(5.42) 



(5.43) 



5.7 Integrable systems related to the g-Krawtchouk polynomials 

The Hamiltonian 



^ (l_g-atai)(i^pgaSao) a>o-a?a,-l (1 + pq^ ^0^0+^1 {I - q-Q-o) 

^ (l+pg2a5ao)(l+p^2aSao+l) +^^5 



(l+pg2a5ao-l)(l+pg2a5ao) +(5-44) 



/ pg^Sao-aJai+l^i ^ p^2aSao+a*ai+l^(i _ ^a^ao+l^^l _ g-ajai^^^ ^ pg^S^o) ^ 



3-0^11 



p^ajao-a'ai+l^]^ + pg2aSao+a5'ai+l')(-]^ _ ^aSao+l)(l _ q-^t^i)(^l + pg^S^o) 



(1 +M2aSao)(l + pg2aSao+l)2(i + pg2a*ao+2)(-a*ao + l)ajai 

where < g < 1 and p> 0, has the spectrum 

a(H,) = {q-':i = 0,l,...}. 

Eigenspaces G '''(H/) are 

= span||£';,jv) : TV = Z + n, n = 0, 1, . . . |, 

where for Ei^n = q~^ 

N 



\Ei,n) = J2^n{Ei,N;p,N;q)\n,N-n) 



(5.45) 



(5.46) 



(5.47) 



n=0 



and 

Kn{Ei^N;p,N;q) 



{-p, q-^; q)n{l + pq''^)p^ q'^''^') 



{q, -pg^+i; g)„(l +p){-pq; q)N{-pq-^Yq'^ 



2 3<P2 



q '-,q ',-pq" 
9-^,0 



a; g (5.48) 
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r 1 ^ 

are the q'-Krawtchouk polynomials. The weight function for the orthonormal system < Kn(Ei f^;p,N;q) > 

I ' i n=0 



IS 



1 



{Ei^n\Eln) {q]q)i 

We can summarize that the spectral decomposition of the interaction Hamiltonian is 

oc N 



JV=0 1=0 



-I \Ei,n) {Ei,n\ 
{Ei^n\Ei^n) 



(5.49) 



(5.50) 



5.8 Integrable systems related to the g-Hahn polynomials 

The Hamiltonian 



Hr = 1 



(1 - a/3g2a5ao)(l _ Q,^^2aSao+l) 



(1 - a/3g2aSao+2)(-;L _ ^^^2a-ao + l) 



+ 



+ 



Q,^-aJai+l(^ - ga5ao + l)(l _ Q,^^2a*ao+a5;ai+2) _ ^^a'ao + l) 

' (1 - a/3g2aS'^o+i)(i _ Q,/3g2aSao+2)2(i _ a/3,?2ar";ao+3^ (-a*ao + l)a^ai 



(5.51) 



'(1 - g-'»^'^i)(l - ag^Sao+i)(i _ a/3g'^5ao+i)aoai + 



Q,^-a*ai + l^^ - ga*ao + l)(l _ Q,^^2a* ao+a^ aj +2) _ ^^a*ao + l) 



(1 - a/3g2aSao+i)(i _ ^^^2aSao+2)2(i _ a/Jg^aSao+s^^agao + l)a*ai 



X ^(1 - q-<^^){l - ag^Sa°+i)(l - a/3g^oa°+i), 
where 0<g<l, 0<Q!< q~^ and < /? < q~^ has the spectrum 

a(H,) = {g-':i = 0,l,...}. 

Eigenspaces H}, g^' e (t(H/) are 

= span||i;i,jv) : A/' = Z + n, n = 0, 1, . . . | 

where -E;_jv = 



AT 



\Ei,n) = Yl Qn{Ei,N; a, p, N\q) \n, N - n) 



n=0 



and 

Qn{Ei,N;a,P,N\q) = 



^ {(3q-q)N{aq)'' g-^; g)„(l - a/?g2»+i)gA-»-fe) 

V {apq'\q)N {q,al3q^+\ Pq-q)n{l ~ al3q){-aqY 

are the g-Hahn polynomials. The appropriate weight function is 



3f2 



aq,q 



-N 



1 



{aq,( 



-N. , 



.g ,q)i 

{Ei,n\Ei^n) ~ {q,p-^q-^:q)i 



{apqr . 



We can summarize that the spectral decomposition of the interaction Hamiltonian is 

oo N 



JV=0 1=0 



-I \Ei,n) {Ei,n\ 
{Ei^n\Ei^n) 



(5.52) 



(5.53) 



(5.54) 



(5.55) 



(5.56) 



(5.57) 
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5.9 Integrable systems related to the Dual g— Krawtchouk polynomials 

The Hamiltonian 

Hj = (1 + c)g-"*"i + 



y (a^ao + Ija^ai 

a.Qa.1 



(a*ao + l)ajai 

where 0<5<l,c<0, has the spectrum 

(t(H/) = {q-^ + cq'-^ ■.l^O,l,...N, TV = 0, 1, . . .}. (5.59) 

The eigenvalues i?;,jv — 9 ' + c<i^^ , Z = 0, 1, . . . , TV of the reduced Hamiltonian Hjv depend on N , 
and therefore the eigenspaces of H/ are one dimensional given by the vectors 

N 

l^^i,Jv) =^/^n(-Bi,iv;c,iV|q) \n,N~n), (5.60) 

Tl = 

where the coefficients 



K4E,„C,N\q) - 302 ^-.^0 



9;<zj (5.61) 

r 1 ^ 

are the dual g-Krawtchouk polynomials. For each fixed N the finite family < Kn{Ei,]\[; c, N\q) > 



N 

is 

•1=0 

an orthonormal system with respect to the weight function 



1 _(c,-,,-;,).(l-c,--)^_,^,,,_,^ ^^^^^^ 



{Ei^n\Ei,n) {q,cq;q)i{l-cq-^) 
We can summarize that the spectral decomposition of the interaction Hamiltonian is 
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